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Abstract 

We present a calculation of the renormalized heavy-light and heavy- 
heavy currents in HQET at order ©(I/ttIq). 

1 Introduction 

In im we have presented a new calculation of the 0{a) renormalized La- 
grangian of Heavy Quark Effective Field Theory (HQET) at 0{l/mQ). This 
result has been confirmed using different calculational strategies in p[ (see 
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also 0). In this note we proceed with the results of a calculation of the 
renormalized heavy-light and heavy-heavy currents in the same order. The 
goal is a systematic calculation of the matrix element 

{D{v')\c^,{l-^,)b\B{v)) (1) 

up to and including O^l/mg), including leading order QCD radiative cor- 
rections. 

A full calculation of the matrix element (|I]) of the heavy-heavy current 
at non-zero recoil is technically feasible, but still of limited phenomenolog- 
ical value. The appearance of the second independent velocity vector v' 
leads to a proliferation of operators whose matrix element would have to be 
parametrized. This results in a substantial decrease in predictive power. In- 
stead we will follow 0] and concentrate on the zero recoil case v = v', which 
should be considered as the lowest order of an expansion of the matrix ele- 
ment (|I|) at the symmetry point uj = vv' = 1. In this case Luke's celebrated 
theorem guarantees that there are no 0{l/mQ) corrections to the current, 
but there are non-trivial contributions at 0{l/mQ) ^. 

This note is organized as follows: in section ^ we define our operatorbasis 
for the effective lagrangian, heavy-light and heavy-heavy currents up to and 
including 0{l/mQ). Our result for the anomalous dimensions will be pre- 
sented in section followed by a discussion of various consistency checks in 
section |^. In section ^ the numerical solution of the renormalization group 
equation (RGE) for the currents at 0{l/mQ2) is presented and is discussed 
with respect to phenomenological applications. 

2 Operator basis 
2.1 Heavy-light currents 

We start with the current up to C(l/mQ) and the operators in the effective 
lagrangian needed for the construction of time-ordered products in 0{l/mQ). 
The kinetic and chromomagnetic operators in the effective lagrangian of di- 
mensions 5 

C,Wc/b ^ ic/.^^^Yhf, Oj^)^/' = Ihfa^'^G.^hf (2a) 

and 6 

^{2)c/b ^ JfJb^J^f^i^i^l^fYJJt^ffJb^ Q^c/b ^ lc/b^^t.v^j^^^^^^^j^YD^h''J^ (2b) 
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are labeled by a flavor index for future use. We choose the operatorbasis for 
the currents of dimensions 3 and 4 as 



Ohl,i = clv'-K, oP^,2 = CL^K (3a) 



and 0« 



—CLiD^^h^, CLiD^hy, —CL{iv D)v'^hy, —CLtpv^hy, 



(3b) 



In (|3bD the label i runs from left to right. Because of the momentum transfer 
at the weak current, we have to distinguish derivatives D = d + igT°'A°' 
acting to the left from derivatives D = d — igT°'A°' acting to the right in 
the operator basis. It turns out that not all operators that contain two of 
these derivatives are independent. In fact there are arrow identities (AI), 
that we have to take into account: 

—i'^ai^p + i'^ pi'^ a — i*Dai'^/3 + i^ i3i*Da =0 (4a) 

iDctiDp — iDaiDp — iDpiDct + iDfjiDa =0 (4b) 

iDaiDp^iDaiDj^ — iDpiDa — iDpiDa =0 (4c) 

cl iDaiDj3 + iDpiDa + iDaiDp + iDfjiDa hy =idaidf^CLhjj (4d) 

The identities (|^) are operator identities and we stress that their derivation |0] 
is independent from the applicability of partial integration. The identity ( ^d]) 
results in additional contributions to the anomalous dimensions from the di- 
vergent parts of the two point function with dimension 5 operator insertions. 
The AI and the equation of motion (EOM) for the heavy quark {ivD)hy = 

(2) 

reduce the full set of 56 local operators to 26 operators Cji/^j, which we 
choose as 



(i?lD)V, iDiDvf", -iDi'^v'', i'^i'^v'', 

-ivDi'^v^, iv3i'pv'^, ivDi^pv'^, -i*Dai^f3cr 

i'^j'^f^a'^f^v^', -ivDi'^i', ivSi^t", i^D^vD, 

{ivDf-f^', i'Di'DY, -i^Di'^'yi', i'^i'^Y, 

-i'Tpi^t', i'pi'^'', -i'DH^, i'^'^i'p, 

i'D^'i'p, -ivDilSaCT'^'', iv3i^c>(T°^, ivDi^D acr"'' 

-i*Djl^pv^ie''^"'f', il^jl^pv^ie''^"'f', 



(5) 
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where the cl on the left and the hy on the right have been suppressed. 
Since iric ^ 0, the full basis contains 18 operators proportional to as well 
as m^, which can be eliminated using the EOM for the light quark 



rricCR 



(6) 



In |T|] the operator basis could be divided naturally in four classes: vanishing 
by the equations of motion (EOM) or not, local or time-ordered product. 
This time we have to consider another class of unphysical operators, that are 
related to the other operators by contraction identities ( CI). These identities 
arise in time-ordered products with (if Z})-operators acting on internal heavy 
quark propagators 0. They are independent from the external states and 
can be illustrated graphically as follows: 

[. . . A{iv'^)K][KB . . oc T[. ..Afi...] 
=^-=^(i^^)==^H®-^ cx =^-@)--*^ 

Even though their derivation proceeds in terms of unrenormalized operators, 
the CI hold true under renormalization in the MS'-scheme, which we haved 
used. The full basis contains 56 time-ordered products of lower order con- 
tributions to the current with operators from the effective lagrangian of the 
6-quark. After application of the CI and the EOM for the heavy quark, we 
are left with 26 time-ordered products. They are classified as follows: 

• double-insertions from the effective lagrangian in (9(l/2m;,): 



-^HL,ijk 



0(0) 



(1)6 QWb 



t,J,k= 1,2 (7) 



The prefactor accounts for insertions of identical operators. To avoid 
double-counting the additional constraint j < k has to be imposed. 



single-insertions from the effective lagrangian in 0{l/{2mi,) 



q-(m)b 

^HL,ij 



0(0) Q 



(2)6 



1,2 



(8) 



mixed insertions of dimension 4 currents and the effective lagrangian 
in C»(l/2m6): 



^(11)6 
^HL,ij 



(l)b 




(9) 
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2.2 Heavy-heavy currents at zero recoil 

The construction of the operator basis contributing to the heavy-heavy cur- 
rent proceeds in analogy to the heavy-hght current. Since the c-Quark is now 

static, the tensor structures of flavor changing operators are projected from 
the left by P+ and we have to consider additional time-ordered products with 
operators in the effective lagrangian of the static c-Quark. To lowest order 
only two operators contribute: 

KP^^v'K KP^i'K (10) 

At order Oil/rriQ), all operators sandwiched between physical states van- 
ish by Luke's theorem and may be discarded. At order Oil/mq) the full 
basis contains 14 local and 100 nonlocal operators. After application of CI 
and EOM for the heavy quarks, we are left with 46 operators. 

The action of spin and flavor symmetry in the matrix elements of time- 
ordered products reveals 20 additional relations. Using flavor symmetry, 
for example, the matrix elements of T^^lf2 and T^^l^^ between mesonic 
states [D) and \B) are described by the same reduced matrix element mul- 
tiplied by Clebsch-Gordan coefficients 

Tr {a^^McPtv^P>''^Mb] , Tt {a^pMcCr''''P+P+v''Mb} (11) 

with P^v'^, a'^^ and M.c/h representing the Dirac-structurc of the current, 
the chromomagnetic operator and the external states respectively. Since the 
former is sandwiched between projectors P/, it conserves spin symmetry and 
may be moved around freely inside the traces. Consequently the traces are 
equal and so are the matrix elements. Other relations among matrix elements 
can be established analoguously. 

Furthermore, since spin symmetry allows us to rewrite all mixed time- 
ordered products in terms of EOM-operators, they are removed from the 
basis applying the CI. 

The basis is reduced by these relations to 6 local operators 

{iDfv^, {iDf^", ipiD", iD^ip, iDJDfsia'"^v^, iDJDpVjie'"^'^^ 

(12) 

taken between static fields h^P^ and and 20 time-ordered products 

^(011)fc6| ^(Oll)cc ^(011)c6| ^(02)6 ^(02)c o\ 

■^HH,ijk\j<ky ^HH,2j2i HH,ijk\j<k for i=l, ■'■HH,iji -^HH,22j K^'^) 

where all indices run from 1 to 2. 
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3 Anomalous dimensions 



As in HI, we have calculated the anomalous dimensions in the background 
field gauge with gauge parameter ^. For the physical operator bases presented 
here, all ^-dependence has to drop out, which provides a powerful cross-check. 
In the case of heavy-light currents, the anomalous dimensions in leading 
and subleading order are well known P, Il0|, ITT] and are reproduced by our 



calculations. Our result in 0{l/ml) can be written in block form, separating 
local operators from time-ordered products: 



.(2) 

Thl 



^{11)6 



(11) 



(11)6 ^(02)6 



(02) 



HL 



.(11) 
Id 



/n,2 



HL 



(Oll)feb 






-.(02) 

Id 

^(011) 



HL 




.(Oil) / 

Id I 



(14) 



Weinberg's theorem guarantees that local operators need no non-local coun- 
terterms. This accounts for the vanishing of seven submatrices in (|T^. There 
are, however, nonlocal contributions 7i°2^'' ^ind 7^°3^'' that arise from the 
renormalization of the lagrangian or the lower order current. This situation 
is depicted diagramatically in figure 

Our result for the heavy-heavy currents at the same order is of similiar 
structure: 



.(2) 

Thh 



^hh 

^(02) 
■^HH 
^(011) 
■^HH 



( 



0(2) 
^HH 

.(2) 
Td 
^(02) 



^(02) 
■^HH 



.(02) 

Td 

/n,2 



^(011) 
■^HH 

^ 


.(Oil) 

Td \ 



(15) 



The entries for all submatrices 7 are available from the authors on request. A 
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Figure 1: Nonlocal contributions to the renormalization of Tjfl^^^'^ . The 
operators of the effective lagrangian in (9(l/mb) and the current in (9(l/m°) 
are represented by open squares and triangles respectively, whereas filled 
symbols represent the corresponding operators of the next dimension. 



typical example is given by the anomalous dimensions of the local operators: 



.(2) 





/o 














o\ 






















Ca 








-1 


1 








4 








1 


-1 






















-2 







Vo 














-2/ 



(16) 



The observation, that the kinetic and chromomagnetic operators (the upper 
left and the lower right 2 x 2-block) renormalize as their counterparts in the 
effective lagrangian in 0{l/mQ), is easily explained by the fact that both 
operator bases coincide at zero recoil up to parity violating contributions. 



The anomalous dimensions have been calculated with the help of FORM ||T2 
In the background field gauge, the divergent part of the three-point function 
has to be calculated to derive the anomalous dimensions. The latter re- 
ceive additional contributions from the two-point function only in the case 
of heavy-light currents as a result of the identity (^). 
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4 Consistency checks 



In section |] we have used several identities to reduce our operatorbasis. How- 
ever, we have to verify that this reduction is compatible with the RGE. Sup- 
pose that a linear combination 

0^ = Y,^^0^ (17) 

i 

of operators does not contribute to the matrix element under consideration 

{D{v)\O4B{v)) = (18) 
In general, this relation will not be stable under renormalization since 

i i ^ i,j 

and 

{D{v)\Os.\B{v))y^0 (20) 

unless jij satisfies further conditions. 

For a consistent reduction of the operatorbasis we have to start with the 
anomalous dimensions of the full operatorbasis. Whenever an operator iden- 
tity (|17[), e.g. CI or EOM, is applied to reduce the full operatorbasis, we have 
verified its compatibility with the RGE. Also, the relations between matrix 
elements of heavy-heavy currents have been shown stable under renormal- 
ization. 

In the background field gauge, only the anomalous dimensions of opera- 
tors that are manifestly invariant under gauge transformations of the back- 
ground field have to be calculated. The Ward identities of the underlying 
BRST-symmetry are particularly simple and provide a powerful cross-check 
of our results []I| . From the Ward identities for the background field and the 
gauge parameter dependence of the effective action follows, that the coun- 
terterms proportional to Ca have to be transversal and should not depend 
on ^. Our results pass these consistency checks. 

Additional checks are provided by the construction of relations among 
operators. For example, the differentiation of the lowest order currents 

id.CLTf'K = CLi'D.T^'K + clT^'i^ ,K (21) 
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leaves their anomalous dimensions unaffected and relates them to subleading 
operators. The operators on the right hand side must mix under renormal- 
ization in such a way, that their sum can be expressed by the total derivative 
of the current on the left hand side, which renormalizes multiplicatively. Our 
result is consistent with all relations of this type, including generalizations 
to local operators of higher dimension and time-ordered products. 

5 Solving the renormalization group equations 

5.1 Matching 

As in any effective field theory calculation, we sum the radiative corrections 
by matching the operators at the thresholds and by running their coefficients 
between thresholds. In the case at hand, we have two thresholds fi = rrib 
and n = rrtc. It is convenient to express the matching conditions as matri- 
ces 

{c\-^,{l-l^m^^^^ = Y.M\{c\OuL,m),=^, (22) 

i 

and M"^ 

i 

Since there is only one QCD-operator, M.^ reduces to a vector. The Ohl/hhj 
contain all physical local operators and time-ordered products up to and in- 
cluding dimension 5, that contribute to heavy-light and heavy-heavy currents 
respectively. Consequently the matching matrices receive contributions sup- 
pressed by inverse powers of the c- and b-quark masses up to and includ- 
ing 0(l/m^y^). In addition, the matching of HL-operators with a covariant 
derivative acting on c induces a nontrivial mc/mft-dependence in Ai'^. For 
example, we have to lowest order 

5.2 RG-Running 

For the lagrangian at 0{l/mQ), only five operators do not vanish by the 
EOM |l| and an analytical solution of the RGE is straightforward. In con- 
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trast, the operator bases for the heavy-hght and heavy- heavy currents con- 
tain 52 and 26 operators respectively and an analytical solution of the RGE 
is impractical, to say the least 

In each order I/ttiq, the solution of the RGE is expressed in terms of 
evolution matrices 

^hl/hh(/^, mc/,) = exp (-^ In (^^) " I'nljnu) ■ (25) 

In P^^^ = (33 — 2nf)/12, the number of active flavors takes the values = 4 
and Uf = 3 foT heavy-light and heavy- heavy currents respectively. 
Introducing their direct sum 

TZhl = © 4l © Rni (26) 

'R-HH = -RhH © -^HH © -^HH ' (27) 

the combined solution for the Wilson coefficients at the hadronic scale can 
be written in compact form: 

CHH(Ahadr) = 7^HH(Ahadr, m,) ■ ■ 7^//L(m„ TTlb) ■ . (28) 



The lowest order currents are conserved at zero recoil, therefore = 1 



,(0) 

in (0). Furthermore, we can set = because Luke's theorem renders 
the corresponding contribution irrelevant. 

For practical purposes, it is most convenient to sum the power series of 
the exponential in (^) directly. This series converges fast and any desired 
accuracy can be obtained. This procedure is used in the numerical results 
in table |l|, where a factor of l/(2mc)^ must be attached to the higher order 
coefficients. We also present the tree-level values of the coefficients to expose 
the effect of the radiative corrections. 

We stress that the matrices Thl/hh defective, i.e. do not have a com- 



plete set of eigenvalues. Therefore the analysis of [|I^ is not applicable and 
we can not derive general properties of the RG evolution from the spectrum 
of the matrices. 

5.3 Results 

Let us finally discuss phenomenological applications of our result. 
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From the coefficient C^^u 2 of the lowest order symmetry breaking opera- 
tor in the top row of table |l|, one obtains a value of 1.180 for the short- distance 
coefficient 77^ of the axial vector current which is larger than 1. This con- 
tradicts the established result in |[T^ which is smaller than 1. This apparent 
contradiction is resolved if one takes into account the complementarity of the 
calculational approaches. 

Our calculation uses effective field theory to resum corrections as In mb/rric 
from the running between and m^, while integrates out the b- and c- 
quarks simultanuously. The latter approach resums the powers of 2; = rric/mf,. 
After applying radiative corrections, includes a" ln"~^ z, a" In*^ z ■ z and 
aghiz-z^ , while we have concentrated on a" In" z, a" In" z ■ z and a" In" z-z"^. 
Therefore our result for r]^ should not be taken face value. Instead, our 
terms a" In" z ■ with n> \ are corrections to the result of JTOl . In practice 



this means that the term a. In z ■ z^ should be subtracted from our result an 



the remaining correction should be added to the result of [jTO| . This correction 
is 0.002 and the final result for riA remains below 1. 

Up to (9(l/mc/f,) no form factors are needed to parameterize the matrix 
elements in the effective theory. The lowest order is normalized and Luke's 
theorem states the absence of contributions at 0(l/mc/fe). 

In 0(l/m^yj,) the constraints of spin and flavor symmetry are less pow- 
erful. An expansion of the matrixelements in terms of 12 nonperturbative 
formfactors has been performed in on tree-level. 

Our result can be used to extend this analysis to include short-distance 
corrections. However, Lukes's theorem, which has been used extensively 
in 1^, receives radiative corrections at 0(l/m^yjj). Therefore the analysis 
of HI has to be refined to include these corrections. 



6 Conclusions 

We have presented a calculation of renormalized heavy-light and heavy-heavy 
currents at order OiXjniQ). The anomalous dimensions of the relevant op- 
eratorbasis of dimension 5 pass various consistency checks. We have solved 
the renormalization group equation numerically for the Wilson-coefficients 
of heavy-light and heavy-heavy currents, matching both quantities at the 
mc-threshold. We have presented the numerical values of the coefficients of 
the heavy-heavy current at the hadronic scale. 

The two expansion parameters io — X and ^hadr/^Q are numerically of the 
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c 



HH,1 



'^HH,13 



a 



0.000 



= A 



hadr 

-0.005 



C 



HH,2 



^HHfi 



a 



1.000 



= A 



hadr 

1.180 



^(02)6 
^HH,ll 
M02)b 
^HH,12 
^(02)c 
^HH,22 
^(011)6fe 

^(011)66 

^(Oll)fefe 
^HH,222 
^(Oll)cc 
^HH,222 

Mon)cb 

^HH,211 

^HH,212 
^(011)c6 
'-^H_ff,122 



-0.667 
-1.222 
2.111 



-0.534 
-1.702 
1.825 



-1.444 
0.556 



C 



(2) 

HHAA 



M02)b 
^ HH.21 
M02)b 
^HH,22 
^(011)66 

^HH.112 

^HH,122 
^(Oll)cc 
'~'HH,212 
,(011)cfe 
'_f/H,lll 
,(011)cfe 
'_ffH,112 
^(Oll)cfe 
^HH.221 

^HH,222 



-2.596 
1.001 



1.444 


-1 


832 


1.111 


-2 


747 


0.111 





023 


0.000 


-0 


005 


0.000 


-0 


004 


0.000 


-0 


002 


1.000 





800 


0.000 


-0 


001 


0.000 


-0 


001 


0.333 





262 


0.333 





155 



0.000 
0.000 
1.000 

1.111 

0.111 
0.111 
1.000 
0.333 
0.333 
0.000 



0.011 
-0.002 
0.421 
1.302 
0.072 
0.043 
0.543 
0.386 
0.229 
0.000 



Table 1: Wilson coefficients in tree-level and one-loop approximation. A 
factor l/(2mc)^ is to be attached to the coefficients of higher order operators. 
For the quark masses we use the values nib — 4.5GeV and rric — l.SGeV and 
for the initial value of the running coupling constant we use as{Mz) — 0.120. 
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same order (for currently accessible values oiu). Therefore phenomenological 
applications should consider also pieces of 0{{u! — 1)^) as well as 0{{uj — 
^)^hadr/fnQ)- While the former have been available for a long time |]14[, the 
latter have been inferred by reparametrization invariance arguments, which 
should be used with care Jll. A direct calculation will be the subject of future 
investigations. 

A more detailed discussion of the consistency checks provided by BRST 
and reparameterization invariance will be presented elsewhere 0, together 
with technical details of the calculation of renormalized Lagrangian and cur- 
rents. 
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